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ABSTRACT 


HS is still a frequently reported endemic disease, with outbreaks in Indonesia. HS vaccines dis- 
tributed in Indonesia exhibit various limitations. This study computationally evaluated the B-cell 
epitope of the 37-kDa OmpH derived from the amino acid sequence of Pasteurella multocida from 
the NTT and Katha strains and compared the epitopes of the two strains. Amino acid sequenc- 
es were obtained from NCBI and analyzed for multiple sequence alignment, and homology was 
analyzed using the BLASTp program at NCBI. Epitope prediction was performed using the IEDB 
B-cell epitope and ABCPred prediction tools. The VaxiJen v.2 online platform was used for antige- 
nicity analysis, and IEDB was used for epitope conservancy analysis. The results of the homology 
analysis revealed that local NTT isolates had a high (>95%) identity with the Katha strain and 
isolates from China, India, Iran, Japan, and the USA. The epitope predictions from both methods 
were cross-checked, overlapping epitopes were shortlisted, and only five epitopes were predicted. 
Among the five, one epitope, ALEVGLN, appeared to be antigenic to both NTT and Katha strains. 
The antigenic sequence of 37 kDa OmpH can be used for peptide-based vaccine development and 
immunotherapeutic design. 
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Introduction 


H also known as snoring disease in Indo- 
nesia, is caused by Pasteurella multocida, a 
gram-negative bacterium that commonly infects cat- 
tle and buffalo [1]. 

HS has a wide distribution, especially in tropical 
regions, such as Middle-East, Central Africa, North- 
Fast Africa, South Africa, South Asia, and South-East 
Asia [2]. In Indonesia, 12,058 cases of HS were report- 
ed between July 2007 and December 2019, with cases 
in buffaloes (4,047), cattle (5,809), pigs (2,108), goats/ 
sheep (64), and Equidae (30) [3]. 

The disease causes huge losses due to livestock 
mortality, as well as losses in the meat- and dairy-re- 
lated industries [4]. The economic loss of the livestock 
industry is estimated at 792 million USD per year [5], 
and the mortality rate of HS is up to 100% [6]. Consid- 
ering its socioeconomic impact, HS is classified as one 
of 22 strategic infectious diseases in Indonesia, where 
control measures are coordinated at the central level 
[7]. 

The disease remains endemic in Indonesia, and 
several regions continue to report outbreaks [8, 9, 
10] despite vaccination efforts. Currently, the disease 
is present in multiple areas of Indonesia, including 
Bali, Bengkulu, DKI Jakarta, Java, Aceh, West Nusa 
Tenggara, East Nusa Tenggara, North Sumatra, South 
Kalimantan, South Sumatra, Jambi, Riau, Central Su- 
lawesi, West Sumatra, East Kalimantan, and Central 
Kalimantan [3]. 

The HS vaccines in Indonesia are administered 
once annually in the form of either an oil adjuvants 
vaccine or an aluminum precipitate vaccine, both of 
which are developed using the Katha strain from Bur- 
ma [7]. The most effective option in the market is the 
oil adjuvant vaccine, which can provide immunity for 
up to one year. However, this vaccine has limitations, 
including the high viscosity of the solution, which 
makes injection challenging and can lead to swelling 
and abscesses at the injection site [2, 11]. Furthermore, 
this type of vaccine has a short storage time because it 
is susceptible to damage caused by temperature fluc- 
tuations [12]. 

The OmpH is a prominent protein in the purified 
envelope of P. multocida envelope and was identified 
as an immunodominant porin. Ongoing research has 
explored the feasibility of utilizing OmpH as a subunit 
vaccine in both native and recombinant forms to com- 
bat avian cholera, bovine respiratory illness, and HS 
[4, 13-16]. Vaccination with the recombinant form of 
OmpH, which has a molecular weight of 37 kDa, can 
induce both antibody- and cell-mediated immune re- 
sponses in dairy calves and swamp buffaloes, protect- 
ing HS [4, 17]. Research conducted by Maulana et al. 
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(2018) demonstrated that 37-kDa OmpH from a local 
strain of P. multocida isolated from NTT and Katha 
strains exhibit similar B-cell epitopes [18]. 

Immunoinformatics is a valuable technique for 
identifying new antigenic epitopes that can be used 
to design new vaccines against a variety of infectious 
diseases, including peptide-based vaccines. Given the 
advancement of immunoinformatics technologies, 
more accurate prediction of B-cell epitopes can now 
be made. Compared to strict laboratory studies, em- 
ploying immunoinformatics methods to anticipate 
epitopes and develop peptide-based vaccines mini- 
mizes costs and saves time, whilst also elevating pre- 
cision [19, 20]. 

This study was conducted to computationally eval- 
uate the B-cell epitope of the 37-kDa OmpH derived 
from the amino acid sequence of P. multocida from 
NTT and Katha strains, and compare the epitopes of 
the two strains. The results can be used for further 
evaluation in the development of a peptide-based vac- 
cine candidate. 


Results 
Multiple Sequence Alignment and Homology 
Analysis 


Conserved and varied regions of the 37 kDa 
OmpH from NTT and other sources were evaluated 
through homology analysis and multiple sequence 
alignment as depicted in Figure 1. Specifically, Clustal 
W was employed to assess seven different sequences of 
the OmpH gene from NTT-Indonesia, China, Japan, 
India, Iran, the USA, and Katha strains. Most of the 
alterations were point mutations, whereas there were 
identified insertions at positions 66-71, and deletions 
at positions 60-61, 189, and 196-197. The Entropy 
(Hx) Plot in Figure 2 indicates high entropy, reveal- 
ing a significant variation between each sequence. The 
highest entropy values of 0.95570 were was recorded 
at positions 6. The homology analysis results showed 
that the 37 kDa OmpH asam amino sequence in local 
NTT isolates was significantly similar (>95%) to the 
Katha strain and isolates from China, India, Iran, Ja- 
pan, and the USA (Table 1). 


Prediction of the B-Cell Epitope of the OmpH P. 
multocida isolates of NTT and Vaccines 


Kolaskar and Tangoankar method was employed 
to predict both the NTT isolates and Katha strain, 
possessing 11 potential B-cell epitopes within the 37 
kDa OmpH gene that met the specified threshold val- 
ue. In addition, the ABCpred prediction server was 
utilized in this study for predicting linear B-cell epi- 
tope regions in an antigen sequence through an ar- 
tificial neural network. The results of both methods 
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Figure 1. 
Multiple alignment of amino acid sequences of 37 KDa gene OmpH P. multocida NTT isolates compared with and isolates 
from GenBank 
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Figure 2. 
Entropy (Hx) Plot graphic of multiple sequence alignment of all amino acid sequence of 37 KDa gene 
OmpH P. multocida NTT isolates with vaccines and isolates frim GenBank 


Table 1. 


Homological value of the nucleotide sequences of OmpH P multocida NTT isolate and vaccine 
ae NTT Isolate Katha Strain 
Comparison isolate (Data GenBank) 
QDC35621 (%) QDF60493.1 


Outer membrane protein [Pasteurella multocida]. Accession 


99.30% 98.95% 
Number ABR24803.1 - China 
Adhesive protein, partial [Pasteurella multocida]. Accession 

98.95% 100.00% 
Number ABX58059.1 - India 
Outer membrane protein, partial [Pasteurella multocida]. Ac- 

96.92% 95.24% 
cession Number AAC02257.1 -USA 
OmpH, partial [Pasteurella multocida]. Accession Number 

96.54% 95.53% 
AQM74565.1 - Iran 
Outer membrane protein [Pasteurella multocida]. Accession 

95.85% 94.85% 


Number ABD94067.1 - Japan 


Katha Strain Accession Number QDC35621 98.95% 100% 
2(3)=1532, p =< 0.001 


were cross-checked, and overlapped epitopes were selected and listed in Tables 2 and 3. 
shortlisted. Next, only five epitopes were ultimately 


Computational Evaluation of B Cell Epitope of 37 kDa Outer Maulana et al., IJVST 2024; Vol.16, No.1 
Membrane Protein H (OmpH) Pasteurella multocida DOI:10.22067/ijvst.2024.81545.1237 


IRANIAN JOURNAL OF VETERINARY SCIENCE AND TECHNOLOGY 


Antigenicity Analysis 


Antigenicity analysis was 
conducted using the VaxiJen 
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Table 2. 
Potential shortlisted BCL epitopes of OmpH P multocida NTT isolate using Kolaskar 
and Tongaonkar antigenicity method and ABCpred method 


j Kolaskar and Tongaonkar Amino Acid an@ored Method Start 
test with a threshold value of 0.5 Antigenicity Method Position ? Position 
tables 4 and 5). f the fi 
( b i d 5), Out of thë five DVGVSDYTYFLG 100-111 DVGVSDYTYFLGGINN 98 
epitopes, one appeared to be an- 
tigenic in both NTT and Katha GAYVFSA 136-142 GFTFGGAYVFSADADK 131 
strains with the same peptide GFVVAGL 154-160 RGFVVAGLYNRKMGDV 153 
sequence ALEVGLN and a 
score of 1.3471, SQKYVTVA 177-184 AGYSQKYVTVAKQEKE 174 
: ALEVGLN 223-229 ALEVGLNYDINDKAKV 223 
Conservancy Analysis 
Conservation analysis was 
Table 3. 


performed on 100 sequences re- 
trieved from the NCBI database 
via the BLASTp tool against the 


Potential shortlisted BCL epitopes of OmpH P. multocida Katha strain using Kolaskar 
and Tongaonkar antigenicity method and ABCpred method 


i Kolask T nkar Amino Aci 
NIT isolate. Only ALEVGLN a a me eres ABCpred Method oat 
. . Antigenicity Method Position Position 
was deemed as a possible anti- 
gen when it came to antigenic- DVGVSDYTYFLG 100-111 DVGVSDYTYFLGGINN 100 
ity, while both GFVVAGL and GAYVFSA 138-144 GFTFGGAYVFSADADK 133 
ALEVGLN showed relatively 
i : GFVVAGL 156-162 RGEVVAGLYNRKMGDV 155 
high conservation levels (93%) 
(Table 6). SQKYVTVA 179-185 AGYSQKYVTVAKQEKA 179 
ALEVGLN 223-229 ALEVGLNYDINDKAKV 223 
Table 4. Table 5. 


Prediction of B-cell epitope antigenicity from outer membrane 
protein H (OmpH) amino acid sequences of P. multocida local 
isolates NTT 


Antigenicity test by Vaxijen v 2.0 


Sequences Score* Probability 
DVGVSDYTYFLG 0.2146 Probable Non-Antigen 
GAYVFSA -0.1360 Probable Non-Antigen 
GFVVAGL -0.1978 Probable Non-Antigen 
SQKYVTVA 0.2844 Probable Non-Antigen 
ALEVGLN 1.3471 Probable ANTIGEN 


“threshold score 0.5 


Table 6. 


Prediction of B-cell epitope antigenicity from outer membrane 
protein H (OmpH) amino acid sequences of Katha vaccine strain 


Antigenicity test by Vaxijen v 2.0 


Sequences Score* Probability 
DVGVSDYTYFLG 0.2146 Probable Non -Antigen 
GAYVFSA -0.1360 Probable Non -Antigen 
GFVVAGL -0.1978 Probable Non -Antigen 
SQKYVTVA 0.2844 Probable Non -Antigen 
ALEVGLN 1.3471 Probable ANTIGEN 


*threshold score 0.5 


Conservation analysis of selected peptides from P. multocida 
NTT isolate and Katha strain through the IEDB epitope conser- 


vancy analysis website 
Peptide Sequence 
DVGVSDYTYFLG 
GAYVFSA 
GFV VAGL 
SQKYVTVA 


ALEVGLN 
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%Protein Sequence Identifi- 
cation 


84.00% (84/100) 
74.00% (74/100) 
93.00% (93/100) 
28.00% (28/100) 


93.00% (93/100) 
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Discussion 

The OmpH is an antigenic surface protein in the 
envelope of P. multocida. It has been identified in all 
bovine isolates and is being considered a vaccine can- 
didate by some researchers [2]. Conserved and var- 
ied regions of the 37 kDa OmpH protein from NTT 
and other sources were assessed through multiple se- 
quence alignment. The entropy (Hx) plot graphs were 
utilized to illustrate the variations between sequenc- 
es. The entropy values had a range of 0-1. The more 
varied a sequence, the higher the entropy value. The 
positions with the highest entropy were 61 and 299. 
Sequences resulting in high entropy alignment are 
unlikely to make good vaccine candidates. An entro- 
py value of 0 was observed in regions 1-7, 14-27, 29- 
59, 78-188 191-195, 198-235, 240-248, 254-265, and 
267-298. In the low-entropy alignment, it is crucial to 
consider regions that correspond to very low entropy 
when designing peptide vaccines. This highly con- 
served region across antigens results in the formation 
of antibodies by the peptide, which protects against 
pathogens of different serotypes. Furthermore, this 
region demonstrates several critical functions due to 
its high level of conservation, making it resistant to 
mutation under positive selection pressure. In addi- 
tion, the binding of antibodies to this region deacti- 
vates essential functions [21]. The 37 kDa amino acid 
sequence of OmpH local isolates in NTT showed high 
identity (>95%) with the Katha strain and isolates 
from China, India, Iran, Japan, and the USA, suggest- 
ing that the nucleotide sequences of NTT, vaccines, 
and other reference isolates may share a common an- 
cestor [22]. 

A B-cell epitope refers to a defined region on the 
surface of an antigen that binds to an antibody. There 
are two categories of epitopes, comprising continu- 
ous (linear or sequential) and discontinuous (non- 
linear or conformational) epitopes [23]. Linear B-cell 
epitopes are made up of peptides that can be readily 
used in immunization and antibody production as 
substitutes for antigens [24]. The Kolaskar and Tan- 
goankar methods revealed that both the NTT isolates 
and Katha strain possess 11 potential B-cell epitopes 
in the 37 kDa OmpH gene that surpass the specified 
threshold value in this study. The Kolaskar and Ton- 
gaonkar antigenicity techniques evaluate the antige- 
nicity dependent on the physiochemical properties of 
amino acids as well as the number of experimentally 
confirmed epitopes [25]. 

Another prediction method for B-cell epitopes 
utilized in this study was the ABCpred prediction 
server. This server implements an artificial neural net- 
work to predict the linear regions of B-cell epitopes in 
an antigen sequence. The methodology for this serv- 
er is developed on a recurrent neural network (ma- 
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chine-based technique), which relies on fixed-length 
patterns. The accuracy of this prediction method is 
65.93%. 

ABCpred was utilized to identify potential B-cell 
epitopes from 37 kDa OmpH P. multocida NTT iso- 
late and Katha strain. The results from both methods 
were cross-checked, narrowed down, shortlisted to 
overlapping epitopes, and only five were selected. The 
combination of both methods can enhance the results 
of the prediction accuracy of B-cell epitope [26]. 

The peptides in the epitope prediction were cho- 
sen using the Kolaskar and Tongaonkar antigenicity 
methods and subsequently verified with the ABCpred 
method. They were then tested for antigenicity using 
the VaxiJen test, with a minimum threshold value 
of 0.5. VaxiJen is a pioneering server for the align- 
ment-independent prediction of protective antigens 
[27]. Out of the five epitopes, one appeared to be an- 
tigenic in both NTT and Katha strains with the same 
peptide sequence ALEVGLN, scoring 1.3471. The 
high antigenicity value of this peptide sequence in- 
dicates its potential as an ideal vaccine candidate. A 
higher antigenicity value translates to a greater abili- 
ty to induce the production of specific antibodies by 
B-cells [28]. 

The 37 kDa OmpH peptides with fairly high con- 
servation (93%) comprise GFVVAGL and ALEV- 
GLN, although only ALEVGLN is considered a prob- 
able antigen. In a study conducted by Bui et al. [29], 
conservation refers to the part of a protein sequence 
that contains an epitope that is considered to be at or 
above a certain level of identity. Regions with lower 
variability in the parts of the protein sequence that 
contain epitopes indicate greater epitope uniqueness. 
This indicates the degree of variability or uniqueness 
of the epitope. Therefore, these regions often serve as 
good targets for the development of epitope-based 
vaccines since the targeted epitopes are present in var- 
ious strains of specific pathogens. 


Conclusion 


According to the findings of our bioinformatics 
studies, there was only one epitope of 37 kDa OmpH 
P. multocida both from NTT isolate and Katha strain 
that met the selected criteria of antigenicity and con- 
servancy analysis. The epitope “ALEVGLN” is hy- 
pothesized to be the most potential candidate for 
a seed epitope/peptide vaccine within the 37 kDa 
OmpH. Further studies are needed to investigate the 
conserved regions of OmpH and assess their efficacy 
as vaccines against P. multocida infection. 
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Materials & Methods 
Data Retrieval 


The amino acid sequences of the 37 kDa OmpH P. multocida iso- 
late retrieved from NCBI with accession numbers QDC35621 and 
QDF60493, as well as other reference isolates with accession num- 
bers ABR24803.1, ABX58059.1, AAC02257.1, AQM74565.1, and 
ABD94067.1 were analyzed [31]. 


Multiple Sequence Alignment 


Amino acid sequences obtained from NCBI were analyzed for mul- 
tiple sequence alignment by employing the CrustalW method in the 
BioEdit program. Amino acid sequences of seven different strains 
were included in this analysis to identify the conserved regions of 37 
kDa OmpH [32]. 


Homology Analysis 


Homology analysis was performed using the BLASTp program 
at NCBI [22]. Specifically, the amino acid sequence of the 37 kDa 
OmpH of the NTT local isolate was compared to the Katha strain and 
five other isolates from different countries retrieved from GenBank, 
NCBI. 


Epitope Prediction 

Epitope prediction was performed using the IEDB B-Cell epitope 
prediction tool (http://tools.iedb.org/beell/) with default thresholds. 
The Kolaskar and Tongaonkar antigenicity method was employed, 
which can predict antigenic determinants with approximately 75% 
accuracy. Moreover, ABCpred (http://crdd.osdd.net/raghava/ab- 
cpred/) was utilized, which can achieve an accuracy of up to 65.93% 
using recurrent neural networks with a 0.5 threshold [24]. Predicted 
peptides were acquired and analyzed further for their antigenicity us- 
ing the VaxiJen v 2.0 tools [27]. 


Antigenicity Analysis 

The antigenicity analysis was completed using VaxiJen v.2 online 
platform —_(http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen. 
html). The targeted proteins in VaxiJen2.0 were predicted using the 
auto-cross covariance method with a threshold of 0.5 [30]. 


Conservancy Analysis 


The epitope conservancy analysis tool used in this study was epi- 
tope conservancy analysis (http://tools.iedb.org/conservancy/) at the 
IEDB (Immune Epitope Database) webserver with sequence identity 
threshold set at 100%. Selected peptide sequences from the previous 
step were assessed for their conservation among 100 homologous 
sequences of P. multocida retrieved from the NCBI database where 
BLASTp was performed against the NTT isolate [29]. 
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